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ABSTRACT 

The statistics of catalogued quadruple stars consisting of two binaries (hierarchy 2+2) is stud- 
ied in comparison with triple stars, with respective sample sizes of 8 1 and 724. Seven rep- 
resentative quadruple systems are discussed in greater detail. The main conclusions are: (i) 
Quadruple systems of e Lyr type with similar masses and inner periods are common, in 42% 
of the sample the outer mass ratio is above 0.5 and the inner periods differ by less than 10 
times, (ii) The distributions of the inner periods in triple and quadruple stars are similar and 
bimodal. The inner mass ratios do not correlate with the inner periods, (iii) The statistics of 
outer periods and mass ratios in triples and quadruples are different. The median outer mass 
ratio in triples is 0.39 independently of the outer period, which has a smooth distribution. In 
contrast, the outer periods of 25% quadruples concentrate in the narrow range from lOyr to 
100 yr, the outer mass ratios of these tight quadruples are above 0.6 and their two inner periods 
are similar to each other, (iv) The outer and inner mass ratios in triple and quadruple stars are 
not mutually correlated. In 13% of quadruples both inner mass ratios are above 0.85 (double 
twins), (v) The inner and outer orbital angular momenta and periods in triple and quadru- 
ple systems with inner periods above 30 d show some correlation, the ratio of outer-to-inner 
periods is mostly comprised between 5 and 10 4 . In the systems with small period ratios the 
directions of the orbital spins are correlated, while in the systems with large ratios they are 
not. The properties of multiple stars do not correspond to the products of dynamical decay of 
small clusters, hence the N-body dynamics is not the dominant process of their formation. On 
the other hand, rotationally-driven (cascade) fragmentation possibly followed by migration of 
inner and/or outer orbits to shorter periods is a promising scenario to explain the origin of 
triple and quadruple stars. 
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1 INTRODUCTION 

Formation of binary and multiple stars is a subject of active re- 
search and debate, still rema ining one of the major unsolved issues 
(Zinnecker & Mathieu 2001). Reliable data on the statistics of bi- 
nary and multiple stars in different environments are essential for 
comparison with theory and further progress. So far, binary stars 
have received some attention, but the statistics of higher multiplici- 
ties is still uncertain. The aim of this paper is to review the statistics 
of multiple stars, focusing on quadruple systems, and to relate it to 
the formation scenarios, whenever possible. 

Multiple stellar systems are often considered individually, 
hence the formation of each object can be explained by a unique 
chain of events. Yet, multiples are not rare "freaks" in the stellar 
world. At least 8% of solar-type stars have three or more compo- 
nents (Tokovinin 2008). The nearest star, a Cen, is triple. Multi- 
ple systems are normal products of star formation, and the star- 
formation theory will eventually be able to model their statistical 
properties. 
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A textbook example of a quadruple system is e Lyr - two wide 
visual binaries on a still wider common orbit. Two striking proper- 
ties of this system are the similarity of the masses of all four compo- 
nents and comparable periods of the inner sub-systems, suggesting 
that the formation process was not completely random. We show 
that such quadruples with similar components and inner periods 
are indeed typical. 

One reason why the statistics of multiple stars remains murky 
is the difficulty of accounting for the observational selection. Cat- 
alogued multiples result from random discoveries using a wide va- 
riety of techniques, so quantifying the observational bias seems 
hopeless. A better approach would be to study volume-limited sam- 
ples, but the number of multiples even among nearby G-dwarfs is 
small, while the incompleteness of higher-order hierarchies is still 
significant. So, at present, the only way to look at the statist ics of 
N > 2 stellar systems is by using catalogues. iFekell Jl98lh pio- 
neered the statistical approach to triple stars. The number of known 
multiples has significantly increased since then, making it possible 
to take the first look at systems with more than three stars. About 
half of the known quadruples were discovered during the last 20 
years as a result of systematic radial velocity programs and high- 
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resolution imaging, the latest ones dShkolniketal]|2008h are not 

included in the present study. 

Here we use the updated Multiple Star Catalogue dTokovininl 

1997, hereafter MSC) to study the statistics of quadruple stars com- 
posed of two close pairs orbiting each other. The problem of the 
observational selection is partially addressed by comparing with 
the triple systems, assuming that the selection effects in both cases 
are similar. Although the exact form of this selection filter remains 
undetermined, we can make reasonable assumptions about its be- 
haviour to figure out which statistical features are real. A catalogue 
such as MSC gives a distorted, but still useful information on the 
true underlying statistics of multiple stars. 

We begin with a short review of multiple-star formation pro- 
cesses and their predictions for the statistics (Sect. O. Some pro- 
totypical multiple stars are presented in Sect. [3] to introduce rele- 
vant terminology and to give a feeling of the data. The samples of 
quadruple and triple stars are described in Sect.|4]together with the 
selection effects. The comparative statistics is presented in Sect. [5] 
Section [6] summarises our findings and discusses the implications 
for multiple-star formation. 



2 FORMATION PROCESSES 

A short review of multiple-star formation scenarios is given to put 
the statistics in proper context. The sequence of events, from the 
fragmentation of pre-stellar cores to the N-body interaction and or- 
bit migration, seems well established. Still, theoretical predictions 
vary widely, depending on the initial conditions, physical processes 
involved or modelled, and computational details. This is a rapidly 
evolving field. 

Fragmentation is believed to be the d ominant mech- 
anism forming binary and m ultiple stars {Bonnelll | 200l |; 



Bate. Bonnell. & Bromml 120021: IDelgado-Donate et alj 120041: 
Goodwin, Whitworth, & Ward-T hompsonl 120041 : iMachida et al.1 
2008). [BodenheimeJ d 19781) envisioned a cascade fragmentation 



where a rotating core first collapses into a ring or disk supported 
centrifugally. The ring fragments, and, if the angular momen- 
tum of the fragments is high, they undergo further rotational 
fragmentation, forming a 2+2 hierarchical quadruple. 

Real pre-stellar cores are in the state of complex "turbulent" 
motions, so the cascade fragmentation gives an over-simplified 
picture which, nevertheless, captures the essential physics. A col- 
lapse of an isolated turbulent core of gas initially creates filaments, 
then clumps in the filament (or inter- sections of several filaments) 
form stars or sub-systems (e.g. Krumholz. Klein. & McKee]|2007l : 



iGoodwinT Whitworth, & Ward-Thompsorj|2 004). As an alternative 
to the rotationally-driven fragmentation, prompt fragmentation is 
induced by an external shock wave or by an m = 2 gravitational 
perturbation from another companion. The Larson's velocity-size 
relation V oc R ' 5 means that in pre-stellar cores most of the 
specific angular momentum j is contained in large-scale motions, 
j oc RV oc R . If the core is isolated, this momentum will be 
conserved in the orbital motion of a wide binary or an outer sub- 
system. 

Simulations show that accretion disks can also produce binary 
companions b y fragmentation and can even form multiple systems 
(cf. F ig. 4 in iBonnelll 1200 ll) . IStamatellos. Hubber. & Whitworthl 
1 20071) studied disk fragmentation into several low-mass compo- 
nents, some of which are subsequently ejected. Unlike the first gen- 
eration of stars produced directly by fragmentation, disk-generated 
companions form over an extended period of time, during the 



whole episode of collapse and accretion. IDelgado-Donate et al.1 
d2004l) show, however, that disk-fragmentation companions have 
less chance to accrete as much matter as the first stars and are often 
ejected during subsequent dynamical evolution. 

Accretion. Binary systems produced by fragmentation at 
isothermal stage of collapse are relatively wide, from 10 2 to 
10 4 A U (cf. the scaling formula in lSterzik. Durisen. & ZinneckeJ 
2003). The fragments continue to accrete surrounding gas, so the 
final paramete rs of a binar y (or a sub-system) are determined by 
the accretion. iBatd «|2000) studied the evolution of the separa- 
tion and mass ratio of an accreting proto-binary under the as- 
sumption that there is no angular momentum transport in the ac- 
creted gas. The results depend on the density and rotation profiles 
of the cloud. Typically (but not always) the mass ratio of an ac- 
creting binary increases and tends to unity, while the orbital pe- 
riod decreases. This trend is confirmed by numerous hydrodynam- 
ical simulations. Binari e s with nearly identical compon ents, twins 
dHalbwachs etai]|2003l : lLucy||2006l : IS6derhielmll2007l) . are most 
naturally explained by accretion. Fragmentation occurs also at adi- 
abatic st age of the collapse a nd at higher densities, creating close 
binaries dMachida et al.l2008l) . However, in this case the fragments 
have very low masses. If these fragments do not merge, the param- 
eters of the resulting close binaries are still settled by the accretion. 

N-body dynamics plays a role in a non-hierarchical clus- 
ter of proto-stars that is formed when several fragments fall to 
the common mass centre of a pre-stellar core. The cluster disin- 
tegrates, ejecting some stars and leaving a stable binary or mul- 
tiple system, usually with the most massive stars paired together. 
As the free-fall time of stars and gas to the centre is the same, 
accretion and dynamics pr oceed simultaneously. Pu re N-body dy- 
namics of small clusters dSterzik & Duriser] 1 19981) produces un- 
realistic multiples with small period ratios and la rge outer eccen- 
tricities (cf. Fig. 6 in ISterzik & Tokovininl l2002h . The presence 
of gas does change the character of the process compared to the 
point-mass dynamics. Extensive simulations of fragmenting cores 
including accretion and dyn amical evolution were performed by 
IDelgado-Donate et al.l d2004l) . They find that the chaotic dynam- 
ics disrupts mostly wide and/or low-mass companions with low 
binding energies, but has little effect on the inner sub-systems with 
massive components formed by fragmentation and accretion. Com- 
ponent masses in these surviving multiples are comparable, as in 
the real systems (see below). In the 10 simulated runs, three 2+2 
quadruples were formed within 0.5 Myr from the start of the col- 
lapse, two of them with additional outer components. All those 
quadruples survived subsequent dynamical evolution of the mini- 
clusters. 

Many authors (e.g. ISterzik. Durisen. & ZinneckeJ 120031 : 
iGoodwin. Whitworth. & Ward-Thompsonl |2004|) believe that N- 
body interaction is a necessary step in the formation of multiple 
stars because it helps to reduce orbital separations by an order of 
magnitude compared to the fragmentation scale. We question this 
assertion. As large-scale motions of the core contain significant an- 
gular momentum, the fragments do not necessarily fall onto the 
centre and interact dynamically, but rather may end up in a stable 
wide binary, possibly containing inner sub-systems. This scenario 
is conceptually close to the cascade fragmentation and it can pro- 
duce a hierarchical 2+2 quadruple s ystem like e Lyr with simila r 
components' masses and inner orbits (Delg ado-Donate et al. 120041) . 
The orbital scales of the e Lyr system roughly match the initial 
fragmentation scales, so no further orbit evolution is required. The 
eccentricity of the outer orbit in this scenario should be moderate, 
the multiple system is hierarchical and dynamically stable. 
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Figure 1. Hierarchy of the selected quadruple stars. The periods of sub- 
systems are shown by the horizontal lines, the size of the circles depicting 
stellar components is roughly proportional to their masses. Note that the 
ratios of outer to inner periods are typically large. 

Orbit migration is an accepted theory for "hot Jupiters", 
where the orbit of a planet shrinks by interacting with the gaseous 
disk. It is necessary to evoke some kind of migration for explain- 
ing the origin of close binaries, unless they are formed directly 
by fragmentation at high densities. The eccentricity distributions 
of exo-planets and spectroscop ic binaries are strikingly similar 
(Ribas & Miralda-Escude 2007), suggesting that interaction with 
disks may be important in both cases. Unlike planets, stellar com- 
ponents cannot migrate in low-mass debris disks, but massive ac- 
cretion disks at earlier stages could act in a similar way (type III mi- 
gration, Peplinski, Artvmowicz, & Mellema 2008). Accretion onto 
a binary also can redu ce its perio d, even without any angular mo- 
mentum transport (e.g. lBate1l200fJ) . but it is likely that in this case 
massive disks are formed as well and affect the orbit of the binary. 

There may be several mechanisms of orbit migration. Their 
common characteristic is that the orbital angular momentum is 
transferred to some other body, while the potential energy re- 
leased from the component's approach is dissipated. Apart from 
the disk migration, can didate processes are interactions with a jet 
teeipurth & Aspinll2004l) or with a ma gnetic wind. Kozai cycles 
with tidal friction (KCTF,|Eggleton 2006) can be viewed as yet an- 
other flavor of migration where the angular momentum is deposited 
into a tertiary component while the orbital energy is dissipated by 
tides. 

Statistics of close binaries help to quantify the migrati on ob- 
servationally. It has been established l lTokovinin et al J 120061) that 
not all low-mass spectroscopic binaries (SBs) with periods below 
30 d have tertiary components, so these SBs are produced by a mi- 
gration mechanism distinct from the KCTF. The SBs with tertiary 
companions have shorter periods than pure binaries, but the mass 
ratio distributions of these two sub-populations are identical, hence 
the dominant migration process does not change the mass ratio. 
Many SBs have tertiary components too distant to cause KCTF mi- 
gration, yet their periods are, statistically, shorter than the periods 
of pure SBs, so the migration is likely enhanced by the presence of 
even a distant companion. 

The properties of multiple stars are reviewed with the aim to 
establish which formation processes are dominant. 

3 TYPICAL QUADRUPLE SYSTEMS 

The data for this study are taken from the current version of the 
Multiple Star Catalogue (MSC) jTokovininl [19971) . We begin by 
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presenting seven typical quadruple systems to illustrate the nature 
of the input data and its reliability. FigureQ]depicts the hierarchy of 
these systems, where a hierarchical multiple system is represented 
as a combination of elementary binaries composed either of sin- 
gle stars or of closer pairs. Relevant parameters of the systems are 
gathered in Table [T] with notations explained in Sect. 14. 11 All sys- 
tems are identifie d by WDS codes based on the 2000 coordinates 
dMason et al]|200lh . with additional identifiers provided in the 2- 
nd column. We focus only on quadruple systems composed of two 
close pairs on a common wide orbit (hierarchy 2+2) and leave aside 
quadruples of hierarchy (2+l)+l, i.e. a hierarchical triple orbited by 
a distant companion. Some triples and quadruples considered here 
contain additional components, as illustrated by Fig. Q] The terms 
triple and quadruple generally mean systems containing exactly 3 
or 4 stars. In this paper we use them to denote objects from our 
samples, so some systems called here "triple" contain more than 3 
stars. 

Epsilon Lyrae is a visual quadruple composed of four A-type 
Main Sequence (MS) stars arranged in two pairs AB and CD. The 
component C has been resolved twice with speckle interferometry 
into a 0.19" binary Cab=CHARA 77, but the existence of this sub- 
system remains doubtful because of many negative observations. 
So, we ignore this uncertain couple Cab. However, some compo- 
nents believed now to be single stars may be resolved in the future. 

Alpha Geminorum is a classical multiple star where both 
components A and B of a 445-yr visual binary are close spectro- 
scopic pairs with A-type primaries and low mass ratios. The eclips- 
ing pair of red dwarfs YY Gem (component C) at angular distance 
72.5" (projected distance 1145AU) also belongs to the system, 
making it sextuple. Thus, there are 2 quadruples of (2+2) hierar- 
chy here: one with the outer 445-yr orbit AB and another with the 
outer orbit AB-C and the inner sub-systems Cab and AB. In the 
following, we ignore composite quadruples, like the outer one in 
a Gem, and consider only simple quadruples composed of just 4 
stars (with, possibly, outer components), like the inner system AB. 
In the case of a Gem, we can be sure that there are only 4 stars in 
the AB system, but our sample of simple quadruples can contain 
some systems with yet undiscovered sub-structure. 

88 Tau is a sextuple system with the same structure as a Gem. 
It is cited here to illustrate that such systems are not exceptional. 
The projected distance between the outer components A and B is 
3500 AU, with A being quadruple and B binary. Recently, interfero- 
metric astrometry has permitted to effectively resolve the inner sub- 
systems Aal,Aa2 and Abl,Ab2 an d to determine th eir orientation 
with respect to the 18-yr orbit Aab JLane et al.l2007l) . Interestingly, 
the inclination of the more massive inner binary Aal,Aa2 relative 
to the outer orbit was found to be 143° ± 2.5°, i.e. it is counter- 
rotating. The inclination of the second sub-system Abl,Ab2 is ei- 
ther 82° or 58°, its orbit is less secure. 

41+ 40 Dra is a well-studie d quadruple with F-type MS com- 
ponents (Tokovi nin et alj [20031) . A combination of observational 
constraints makes the existence of any additional companions ex- 
tremely unlikely, so this system is a genuine quadruple. The mass 
ratios of the sub-systems Aab and Bab are nearly the same, 0.93 
and 0.91. The periods of Aab and Bab are different, but the orbit 
of Aab is highly eccentric, so the angular momenta of Aab and 
Bab are similar. The pairs Aab and AB rotate in opposite directions 
as projected onto the sky. Although the outer orbit is still unde- 
termined, we know that it is not co-planar with the orbit Aab, the 
probable relative inclination is either 60° or 160°. 

GJ 225.2 was a visual triple star until recently, when a 
new sub-system CE has been discovered with adaptive optics 
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Table 1. Selected quadruple systems 



WDS(2000) 


Name 




Periods 






Masses, Mq 








Pl 


Psi 


P S 2 


Mi,i 


Mi i2 


M 2 ,i 


M 2 ,i 


18443+3938 


e Lyr 


340 ky 


1804 y 


585 y 


2.31 a 


1.62 a 


1.86 a 


1.70 a 


07346+3153 


a Gem 


445 y 


9.2 d 


3.0 d 


2.76 a 


0.47 m 


2.98 a 


0.24 m 






14 ky 


445 y 


0.8 d 


3.22 s 


3.23 s 


0.59* 


0.58 * 


04356+1010 


88 Tau 


18 y 


3.6 d 


7.9 d 


2.06 * 


1.36* 


1.07 * 


1.06* 






80 ky 


18y 


3.7 y 


3.42 s 


2.13 s 


1.18a 


0.15 m 


18002+8000 


41+40 Dra 


10.3 ky 


1247 d 


10.5 d 


1.39 * 


1.30* 


1.32 * 


1.20* 


06003-3103 


GJ 225.2 


460 y 


67.7 y 


23.7 y 


0.67 * 


0.52* 


0.69 * 


0.20* 


11221-2447 


HD 98800 


345 y 


262 d 


315 d 


0.80* 


0.67 q 


0.90* 


0.19 m 


04325+1731 


GGTau 


45 ky 


200 y 


2700 y 


0.78 * 


0.68 * 


0.11 * 


0.04* 



Table 2. Sample of triple systems (fragment) 



WDS(2000) 


log P L 


logPs 


Mi 


M 2 


M 3 


Plx 




day 


day 


Mq 


Mq 


Mq 


mas 


00046-4044 


6.33 


3.10 


0.06 a 


0.06 * 


0.26 * 


76.9 


00057+4548 


7.85 


5.74 


0.65 a 


0.47 a 


0.37 a 


85.1 


00063+5826 


4.59 


1.68 


0.86 * 


0.31 * 


0.98 * 


49.3 


00125+1434 


3.19 


0.27 


0.78 m 


0.54 a 


0.52 q 


24.7 


00150+0849 


6.73 


-0.08 


1.24 a 


1.15 * 


1.13 * 


12.5 



(Tokovinin et al. 2005). The period of CE was estimated from the 
astrometric perturbations. So far, there was no radial velocity (RV) 
monitoring of these stars, so additional close components may still 
hide inside. The system is apparently old, as it has a high velocity 
relative to the Sun. A premature orbit of the outer system AB-C 
with Pl = 390 yr was published, but all we can say now is that the 
outer period is of the order of 500 yr and that all three orbits could 
be close to co-planarity. 

HD 98800 is a young, pre-Main Sequence (PMS) quadru- 
ple consisting of two spectroscopic binaries Aab and Bab at 2" 
from each other. The sub-system Bab was resolved by interferom- 
etry, pr oviding its "visual" orbit and a reliable measurement of the 
masses dBoden et alJl2005h . The outer orbit of AB is only prelimi- 
nary, but there are indications that the AB and Bab orbits have only 
moderate mutual inclination. The mass of Ab is still a lower limit 
from the single-lined spectroscopic orbit. The Bab sub-system is 
more massive than Aab and is surrounded by a massive disk caus- 
ing extinction and flux variability, whereas Aab has no such disk. 

GG Tau is another PMS quadruple only 1-2 Myr old. All four 
components are resolved, the masses are estimated from the evolu- 
tionary tracks, with Bb possibly being sub-stellar. The sub-system 
Aab is surrou nded by a ring-like disk , there is a circum-quadruple 
disk as well. iBeust & Dutrevl d2006t> used the disk parameters to 
constrain the dynamics of this quadruple. They conclude that the 
relative inclinations between Aab, its disk, and the outer orbit are 
likely moderate, so the system is more-or-less co-planar. Interest- 
ingly, the pair Bab is close to the limit of the dynamical stability, 
and, if stable, it is probably counter-rotating. GG Tau may disinte- 
grate during further dynamical evolution, ejecting Bb and leaving 
only a triple system. 



4 STATISTICAL SAMPLES 
4.1 Parameters and notation 

Hierarchical multiple systems are described by a large number of 
parameters - component's masses and orbital elements of each sub- 
system. We concentrate on masses and periods as most relevant to 
the formation processes and touch on other characteristics (eccen- 
tricity, mutual orbit inclination) only briefly. The semi-major axis 
is related to the period by the third Kepler's law and is statistically 
equivalent to period, except for a weak dependence on the total 
mass. 

Consistent notation of the parameters is required to avoid con- 
fusion. For triple stars, the outer (long) and inner (short) periods 
are denoted as Pl and Ps, the masses of the inner binary compan- 
ions as Mi and Mi and the mass of the distant tertiary as M3. The 
inner and outer mass ratios qs and qLz are defined as 

M 2 M 3 
gs= M? gL3 = WTM- 2 > (1) 

although an alternative definition q' L3 — M3/M1 is also possible. 
The outer mass ratio can take values gz,3 > 1 if the tertiary is more 
massive than the inner binary (this happens in 5% of the triples). 

The notation for quadruple stars is as follows: the period of the 
outer system is Pl, Psi and Ps 2 are the periods of the inner sub- 
systems, with index 1 referring to the sub-system with the largest 
total mass. The masses of all 4 components are designated as Mi,i, 
Mi j2 , M 2 ,i, andM 2 ,i. 

The mass ratios of quadruple stars can be defined in various 
ways. We use the following definitions: 

Mi, 2 M 2 , 2 M 2 ,i + M 2 , 2 

qS1= MTi' qS2 = M2~i' to4= M M +Mi, 2 - (2) 
One can argue that the outer mass ratios for triple and quadru- 
ple stars are not directly comparable because the latter compares 
two stellar pairs. An alternative definition of the outer mass ratio in 
quadruple stars can be q' L4 , 

Mia 



Ila = 



9L4 



(3) 



Mi,i + Mi, 2 1 + q S 2 

This parameter can be compared with more directly. For ex- 
ample, when all components have equal masses, we have qn = 1, 
but q' Li = qL3 = 0.5. The choice of the most relevant mass ratio 
definition depends on the implied formation process, e.g. qLz and 
qL4 can characterise the fragmentation of the outer sub-system. 

The data on the periods and masses are taken from the MSG 
The periods are either known more or less exactly from the spec- 
troscopic or visual orbits or, otherwise, estimated roughly from the 
projected separations in the case of wide binaries. The components' 



Statistics of triple and quadruple stars 5 



Table 3. Sample of quadruple systems (fragment) 



WDS(2000) 


log P L 


logPsi 


log Ps2 


M M 


Mi, 2 


M 2 ,i 


M 2 ,2 


Plx 




day 


day 


day 


Mq 


M Q 


Mq 


Mq 


mas 


00134+2659 


7.20 


5.13 


0.78 


2.50 a 


1.54 a 


0.92 * 


0.88 * 


8.1 


00247-2709 


4.16 


2.30 


3.41 


0.10* 


0.08 * 


0.08 * 


0.08 * 


132.8 


00316-6258 


6.68 


5.26 


4.21 


3.84 a 


0.40 v 


2.76 a 


1.94 a 


23.4 


00345-0433 


7.18 


5.56 


1.91 


1.70* 


1.40* 


1.14* 


0.10 m 


10.0 


00364-4908 


8.57 


0.57 


1.34 


0.98 a 


0.83 a 


0.78 a 


0.31 m 


22.8 



Notes: 

00134+2659 = HR 40. VB (368y G0III) and SB (6d G7V) at 18". 

00247-2709 = GJ 2005, visual quadruple composed of M5 dwarfs, 0.27" and 0.05" at 1.07" from each other. 

Short estimated periods, but no orbits yet. 
00316-6258 = HR 126+127+125, Beta Tuc, 6 components. The inner quadruple contains 2 VBs (2.4" B9V and 44.7y A2V) 

at 27" . The distant pair at 544" is itself a 0. 1 " VB with a short, but yet unknown period. 
00345-0433 = HIP 2713. VB (1.9" G8II) and SB1 (81d GOV) at 19.5". 
00364-4908 = HIP 2888 = GJ 24AB. Two SBs (3.7d G3V and 22d K0V) at 329", CPM. 



masses are estimated by various techniques, as reflected by the 1- 
character codes. The most common mass codes are a (mass from 
the spectral type), * (mass estimates from orbits, isochrone fitting, 
etc.), q (secondary of an SB2), m (minimum mass of an SB1 sec- 
ondary), v (estimate from the magnitude difference in a visual bi- 
nary), and : (uncertain mass, e.g. in an interferometric binary with- 
out relative photometry where a magnitude difference of 2 is as- 
sumed). The mass code s means that the "component" consists of 
several stars and the sum of their masses is given. The reader is re- 
ferred to the MSC for further details. Obviously, the mass estimates 
in the MSC are quite crude, so the mass ratios discussed below are 
only indicative of statistical trends, rather than well-measured. 

4.2 Triple and quadruple systems 

The data on triple and quadruple stars are listed in the Tables 2 
and 3, respectively. These tables are available in full electronically, 
their fragments are given here. Each system is identified by its 
WDS (2000) code. Additional identifiers are given in the Notes to 
Table 3. More information on each system can be obtained from 
the on-line MScQ 

We selected all "simple" triple systems composed of an inner 
pair and a distant single tertiary, with possibly some other more 
distant components. The systems with uncertain masses (code :) 
are omitted, leaving TV = 724 triples listed in Table 2. 

The sample of 2+2 quadruples from the MSC contains 135 
entries. We exclude from the statistical analysis composite quadru- 
ples, systems where some periods or masses are not known (e.g. 
astrometric sub-systems without computed orbits), leaving 81 sim- 
ple quadruples. Periods, masses, and mass codes are given in the 
electronic Table 3, together with the notes on individual systems. 
We keep in the sample 3 systems (04226+2538, 11395-6524, 
17146+1423) where one of the 4 masses is poorly determined (code 
:), but exclude some PMS m ultiples with yet unc ertain masses like 
those recently discovered bv lCorreia et al. (2006). Three quadruple 
systems from Table 3 (SZ Cam, n Ori, QZ Car) belong to open 
clusters. 

Both samples are quite heterogeneous with respect to masses, 
ages, and reliability of the data. The most massive is the 40 Mq 
primary component of the quadruple QZ Car (10444—6000), while 

1 http://www.ctio.noao.edu/~atokovin/stars/index.php 
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Figure 2. Top: Cumulative distributions of the primary masses: M\ t \ in 
81 quadruples (full line) and Mi in 724 triples (dashed line). Bottom: The 
number of triple and quadruple systems within given distance d. The dotted 
lines show the d 1 ' 7 trends with coefficients 0.3 and 0.06. 



the red dwarf primary in the quadruple GJ 2005 (00247-2709) 
has the smallest mass of 0.1 Mq. Figure [2] shows the cumulative 
distributions of the primary-component masses in the samples of 
triple and quadruple stars. In both cases about 80% of primaries 
have masses between 0.8 and 5 Mq, with the median mass of 
1.7 Mq. We tried to split the quadruples into low- and high-mass 
parts around this median and found the statistics to be similar. 
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In a reasonably complete nearby sample the number of objects 
within a given distance d would be proportional to d 3 . The number 
of triple and quadruple stars in our samples is rather proportional 
to d 1 ' 7 because low-mass objects are located closer than high-mass 
ones (Fig. [2] bottom). Such weak dependence on d is expected for 
magnitude-limited samples. We deduce from the ratio of the coeffi- 
cients in the d 1 ' 7 fits that triples are 5 times more frequent than 2+2 
quadruples. We recover the N cc d s law out to d < 40 pc if only 
primaries more massive than 1 Mq are retained. The median dis- 
tances to triples and quadruples are lOOpc and 70 pc, respectively. 

4.3 Selection effects 

Multiple stars are discovered by different techniques or their com- 
binations, each with its own biases. Moreover, the choice of the 
surveyed objects is random rather than systematic. As a result, the 
current knowledge of multiplicity is incomplete even in the solar 
neighbourhood. Here we study the statistics of the catalogue, not 
the statistics of real multiple stars in the sky. 

It is standard practise to adopt some models of observational 
selection for deriving true statistics from catalogues or biased sam- 
ples. This approach is justified for some well-defined samples (e.g. 
iTokovinin etal] f2006). but it is problematic for compilative cata- 
logues resulting from random discoveries, such as the MSC. A se- 
lection function f is the ratio of objects in the sample to their true 
number in the same stellar population. Here we intentionally leave 
the selection function undefined and do not attempt to study true 
distributions of periods and mass ratios. However, we investigate 
the correlations between the parameters of the catalogued systems 
and try to figure out the influence of selection effects qualitatively. 
Two reasonable assumptions are made, (i) The discovery of a sub- 
system depends mostly on its period P and mass ratio q, so, to the 
first order, the selection function has the form f(P, q). (ii) The se- 
lection function f(P, q) is smooth, without features favouring or 
disfavouring particular periods or mass ratios. The selection func- 
tion acts as a filter superposed on the observed distribution. If there 
are correlations between the parameters of real multiple stars, they 
are likely to be seen in the catalogue as well. Sharp features in the 
real distributions will be recovered in the catalogue, too. 

The integral of f(P, q) gives the completeness of the sample 
fo. The MSC certainly misses many low-mass companions. How- 
ever, a bold assumption that we see only the "tip of the iceberg" 
and that stars are surrounded by swarms of yet-undiscovered low- 
mass compa nions meets with problems , at least for nearby dwarfs. 
For example. lEggenberger et al] 12007) find that only 15% of stars 
previously considered as single and surveyed with adaptive optics 
have low-mass companions at separations above 7 AU. Precise ra- 
dial v elocities of 93% of "single" dwarfs are stable dNidever et alj 
2002). It is therefore reasonable to assume that the completeness of 
the MSC for nearby dwarfs is of the order fo ~ 0.5 or higher. 

The same techniques are used for discovering inner sub- 
systems in quadruple and triple stars, so the selection functions 
f(Ps,qs) for triples and ,f(Psi,qsi) for quadruples should be 
similar. Unfortunately, we cannot assume that the selection func- 
tions J{Pl,<Ilz) for triples and /(-Pl, <?L4) for quadruples are 
the same. Discovering a faint tertiary can be easier than find- 
ing that this tertiary is in fact double (cf. the system GJ 225.2). 
A systematic study of tertiary components has converted many 
of them into additional pair s, promoting triples to quadruples 
dTokovinin & Smekhovll2002l) . It has been shown that 30% of ter- 
tiary companions believed to be single are in fact binary. 

As mentioned above, only in few multiple systems we can be 



sure that all components are actually known. When a new compo- 
nent is found in the future at some intermediate level (e.g. orbiting 
a close pair with a more distant and already known tertiary), our 
current knowledge of the system's structure is obviously wrong, as 
the values of periods and masses in the Tables 2 and 3 correspond to 
some other hierarchical levels. When a new component is found in 
a close orbit around some star, a simple triple or quadruple is con- 
verted into a complex one. Hopefully, such incompleteness affects 
only a fraction of the systems and the distortion of the statistics 
caused by these missing levels is tolerably small. It can be reduced 
by restricting the samples to nearby, better-studied systems. 



5 COMPARISON BETWEEN TRIPLES AND 
QUADRUPLES 

Periods and masses of the companions are determined by the for- 
mation and early evolution, therefore we study the statistics of these 
parameters. Nevertheless, a triple is characterised by 5 numbers (2 
periods and 3 masses), a quadruple - by 7 numbers. Looking for 
correlations in this multi-parameter space is not easy, so we give 
below only some plots which appear to be most relevant. We begin 
by looking at the periods in the inner and outer sub-systems. Then 
the inner and outer mass ratios are compared. Finally, we explore 
correlations of the inner and outer mass ratios with respective peri- 
ods, the correlations between the orbital angular momenta, and the 
properties of quadruple systems. 

5.1 Outer and inner periods 

Figure [3] plots the triples and quadruples in the period-period dia- 
gram. Throughout the paper, the periods are measured in days and 
plotted on the logarithmic scale. Hier archical system s are dynami- 
cally stable when Pl/Ps > 5 (e.g. Eggleton 2006), this limit is 
shown by the dashed lines. As the periods of wide systems are 
estimated only crudely, some points can fall slightly below this 
line. The distributions of inner and outer periods in 1 85 triple stars 
within 50 pc and in all quadruple stars are plotted in Fig. [4] 

The inner periods in both triples and quadruples are dis- 
tributed in the same way, clustered in two groups: close (log Ps < 
1.5) and wide (log Ps > 4), with a partially filled gap in-between. 
As the discovery techniques for the close sub-systems are predomi- 
nantly spectroscopic and for the wide sub-system mostly visual, the 
gap could be attributed to the selection effect, but it could also be 
real. Note that the inner and outer periods in triple stars show some 
correlation for Ps > 30 d. Mo st points are lo cated in the band 
5 < Pl/Ps < 10 4 (see also lTokovinirill2004l) . while the space 
above this band is almost empty, contributing to the gap in the dis- 
tributions of the inner periods. This correlation must be genuine 
because the inner and outer sub-systems are discovered indepen- 
dently of each other. The correlation does not hold for short inner 
periods Ps < 30 d. 

Tidal interaction in multiple stars (KCTF) shortens the in- 
ner periods, leading to the excess of su b-systems with periods 
of few days (Fabrv ckv & Tremainell2007l) . There is a strong peak 
at logPs = 0.25. ..0.75 in the distribution of inner periods in 
the complete sample of triple stars. However, it can be a selec- 
tion effect produced by the eclipsing systems which are discov- 
ered at long distances, because this peak almost disappears for 
nearby triples within 50 pc (Fig.[4]l. The distribution of the inner 
periods in the homogeneous radial-velocity survey of triple stars 
(ITokovinin & Smekhovll2002T) is nearly flat out to Ps ~ 60 d. So, 
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Figure 3. Comparison of the outer and inner periods in triple (top) and 
quadruple (bottom) stars. All periods are in days. For quadruples, Pgi are 
plotted as diamonds, P§2 as crosses, the points belonging to the same sys- 
tem are connected by dotted lines. The diagonal lines delineate the approx- 
imate dynamical stability limit Pl/ Ps > 5. 



KCTF can produce only a mild depletion at Ps ~ 10 d but cannot 
explain the whole gap. 

Note that the period distribution in the volume-limited sam- 
ples of dwarf binaries is a rising function in the logP inter- 
va l from to 3, with a possible dip near logP ~ 2 (Fig. 10 
in lHalbwachs et ai]|2003l) . The rise is even steeper for non-triple 
(pure) binaries, because all SBs with P < 3 d are triple, while the 
propo rtion of triples decreases at longer periods dTokovinin et al.l 
2006). The obvious conclusion is that the migration mechanism 
produces shorter orbital periods in triple stars than in pure bina- 
ries. It is then plausible that the gap in the inner-period distribution 
in triple and quadruple stars is real and is caused by migration to 
shorter periods, rather than by the observational selection. 

The distributions of the outer periods of quadruples and triples 
are different. All quadruples except two have log Pl > 3.8 (Pl > 
17 y). Moreover, there is a concentration of points near the lower 
limit log Pl ~ 4, not seen in the triple systems where the distribu- 
tion of Pl is smooth and extends to shorter periods. It is shown be- 
low that tight quadruples with log Pl < 4.5 are also distinguished 
by larger outer mass ratios ql4 and more similar inner periods, so 
the reality of this special group of objects is based on more than 
just outer periods. They are further studied in Appendix lAl 




Figure 4. The distributions of the inner (top) and outer (bottom) periods in 
triple and quadruple stars. The bars refer to the period bins of 0.5 dex width 
and are displaced from the bin centers to avoid overlap. In quadruple stars, 
both inner periods are merged. 



The quadruple system with the shortest outer period of 355 d is 
VW LMi, it combines a 0.5-d ay eclipsing pair wit h a 8-day binary, 
all composed of dwarf stars jPribulla et alj|2006l) . The triple star 
with the shortest outer period is A Tau: Pl = 33 d, Ps = 4 d, 
Mi = 7.2M Q . 



5.2 Outer and inner mass ratios 

It appears that there is no correlation between inner and outer mass 
ratios neither in triples nor in quadruples (Fig.O. The distributions 
of qs in triples and qsi , qs2 in quadruples look quite similar. Both 
show some excess of sub-systems with nearly identical m asses, 



-sy st 

twins ( Halb wachs et ai]|2003l ; [Lucvll2006l ; I Soderhi elnj|2007h . 

In triple systems where the mass of the outer companion 
equals the mass of an inner companion, Ms — Mi or M3 = M2, 
there is a relation between inner and outer mass ratios: qL3 = 
1/(1 + qs) and qL3 = 9s/(l + qs), respectively (dotted lines 
in Fig. [5] top). There seems to be a concentration of points along 
these lines, suggesting that the phenomenon of twins may extend to 
the outer sub-systems of triple stars in this strange and un-explained 
way. 

The fraction of triple systems where the outer companion has 
the smallest mass, Ms < M%, is 332/724=46% for the whole sam- 
ple and 71/184=39% for triples within 50pc. The median outer 
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mass ratio is qL3 = 0.39, 81% of triples have qLz > 0.2. 
Thus, tertiary companions tend to have masses comparable with 
the components of the inner binary. Triple systems resulting from 
the N-body decay, on the contrary, have o nly low qL3 (e.g. Fig. 
4 inlPelgado-Do nate. Clarke. & B ate 2003) and do not match real 
triples in this respect. Although the observational selection does 
favor multiples with comparable masses, the difference with the 
simulations is too strong to be explained entire ly by the selection. 
The survey of PMS stars bv lCorreia et al.l ( 120061) confirms that low- 
mass tertiaries are indeed rare, not just missed. 

The distributions of the outer mass ratios in quadruples and 
triples, qn and qL3, are different. A strong tendency to large qn 
is obvious, 97% of quadruples having qn > 0.2. The concentra- 
tion of points in the upper right corner of Fig.|5]corresponds to the 
quadruples where all component masses are similar. Simulations 
of multiple stars with components selected independently are pre- 
sented in Appendix IB"1 to show how different they are from the real 
systems. 

5.3 Inner mass ratios vs. inner periods 

Figure |6]presents a comparison of mass ratios with periods for the 
inner sub-systems. The statistics of triple and quadruple stars ap- 
pear to be similar. We note the relative paucity of inner periods 



between 30 and 10 4 days (the period gap, see above) and the con- 
centration of points towards qs ~ 1 (twins) on both sides of the 
gap. In both samples the sub-systems with log Ps > 4 tend to have 
more similar masses, qs > 0.5. This trend may be caused by the 
selection (sub-systems with low q are not discovered by visual tech- 
niques), as illustrated by the GJ 225.2 system. We checked that the 
distributions of the inner and outer mass ratios qs and q' L3 in nearby 
triple stars with inner (resp. outer) periods in the same range, be- 
tween 10 4 and 10 6 days, are statistically indistinguishable. On the 
other hand, sub-systems with short Ps and low qs could be formed 
by an alternative mechanism such as disk fragmentation. 

5.4 Outer mass ratios vs. outer periods 

Figure|7]shows the correlation of outer mass ratios with outer peri- 
ods for triples and quadruples. These plots are strikingly different. 
In triple stars, we see no dependence of median qLs (diamonds) on 
the outer perio d, although the dispersion of q^z seems to increase at 
larger periods. Delg ado^Donate. Clarke. & Bate! d2003l) predicted 
that close inner binaries with low-mass tertiaries should be com- 
mon, but there is no correlation between qLa and Ps to support this 
prediction. 

Note the complete absence of quadruples in the lower left cor- 
ner with qn < 0.6 and log Pl < 5.4 and the relative abundance of 
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vs. specific momentum in the inner sub-system(s) js in triple (top) and 
quadruple (bottom) stars. The solid lines denote equality of the momenta, 
the dotted lines show Jl /js = 2.5. 



triples with similar parameters. This difference could be caused by 
the observational selection. Quadruples with qn < 0.6 can escape 
detection because the low-mass pair is much fainter than the main 
pair. One such system with estimated qn — 0.63 is \ Tau B where 
the light of the second binary is not detected at all, wh ile its binarit y 
is inferred from the total mass and infrared colour jTorresll2006h . 
Missed tight quadruples with qn < 0.6 and Pl < 100 y would be 
considered today as triples with qL3 < 0.6. We can establish how 
many of such low-gi,3 triples are actually quadruple by monitoring 
the radial ve locity of faint tertiary components wi t h the detection 
technique of lD'Angelo. van Kerkwiik. & Rucinskil d2006h . 

5.5 Orbital angular momenta 

Orbital angular momentum J of a binary star equals 

J = y/a(l - e 2 )MiM 2 ^G/M, (4) 

where Mj and AI2 are the component's masses, M = Mi + M2 
is the total mass, e is the eccentricity, a is the semi-major axis, 
and G is the gravitational constant. In the following we neglect the 
eccentricity because for a typical e = 0.5 the factor (1 — e 2 ) 1 / 2 = 
0.87 is insignificant in comparison with the large scatter of J (there 
are notable exceptions like 41 Dra). We calculate a from periods 



and masses using the Kepler's law. The specific angular momentum 
j — J/M is estimated for each sub-system. 

Angular momentum is strongly related to orbital period, so the 
plots in Fig.[8]resemble the plots in Fig. [5] The correlation between 
inner and outer periods translates into the correlation between the 
momenta, a tighter correlation for wide multiples. Considering that 
J oc a 1/2 oc P 1/3 , the typical ratio J L /Js ~ 10 translates to 
Pl/Ps ~ 10 3 . 

Relative orientation of the outer and inner orbits gives addi- 
tional clues to the formation and dynamical evolution of multiple 
stars. For triples, it has been established that the angular momentum 
vectors of the inner and outer orbits do show some correlation, the 
average angle betwe en them, being about 50 ° instead of 90° for 
un-correlated spins dSterzik & Tokoviiiii42 002). The correlation of 
orbital spins is stronger at low Pl /Ps ratios (close to the dynami- 
cal stability limit) and disappears at large Pl/Ps- 

There are two simple quadruple stars, 88 Tau and /j, Ori, 
with complete elements of both inner and outer orbits known. In 
88 Tau, the inner orbits are not co-planar to the outer orbit, with 
$ > 90° (counter-rotation) for at least one sub-system dLane et all 
120071) . Similarly, the Aab sub-system in fi Ori is counter-rotating, 
$ = 137° ± 8°, while the orbit of the Bab sub-system i s nearly 
perpendicular to the orbit of AB (Muterspaugh et al. 2008). There 
are other cases where the apparent rotation of the resolved inner 
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and outer sub-systems is opposite, suggesting non-aligned spins. 
On the other hand, some wide quadruples are close to alignment 
(e.g. e Lyr, GG Tau, HD 98800). Although this evidence remains 
circumstantial, it hints on the same trend as in the triples, i.e. the 
alignment of orbital spins for low Pl / Ps and mis-alignment for 
large period ratios. A dedicated interferometric survey of inner or- 
bits is obviously needed to increase the number of multiple systems 
with known sense of relative rotation. 



5.6 Quadruple systems 

In this sub-section, we focus only on the quadruple stars and look 
for correlations between mass ratios and periods. There is no corre- 
lation between the inner mass ratios qsi and qs2 (Fig.[9j. However, 
there is a concentration of points towards qs2 ~ 1, and a less obvi- 
ous preference of sub-systems with qsi ~ 1. There is a cluster of 
1 1 points in the upper right corner, qsi > 0.85, qs2 > 0.85 marked 
by the dotted line in Fig. [9] If the points were distributed uniformly 
in the (qsi,qs2) plane between and 1, we would expect only 
2.2% in this corner, while the actual number is 11/83=13%. The 
fraction of double twins in the simulated quadruple systems with 
independent component pairing is 0.25%. The family of quadruples 
where both sub-systems are twins (double twins like 41+40 Dra) 
thus appears to be distinct from the majority of other quadruples. 
BD— 22° 5866 is yet another low-mass quadruple composed of two 



twins dShkolnik et alj|2008l) . The "definition" of twins q > 0.85 is 
adopted here arbitrarily, guided by Fig. [9] 

Sub-systems with qsi > 0.85 or qs2 > 0.85 are found at all 
inner periods (cf. Fig. [6] bottom). Interestingly, the proportion of 
twins in the sub-systems of low mass seems to be higher than in 
the high-mass ones. Similarly, there are more twins among the sec- 
ondary (less massive) sub-systems. These trends need a confirma- 
tion with larger samples. The origin of twin binaries is not yet un- 
derstood. They could result from the preferential accretion onto the 
secondary that dri ves q toward s one and, at the same time, shortens 
the orbital period (Bate 2000). However, the inner mass ratios do 
not correlate with the inner periods (Fig. [6). Figure|9]demonstrates 
the lack of correlation between the periods and mass ratios in the 
two inner sub-systems belonging to the same quadruple. Therefore, 
the inner periods and mass ratios are determined by different pro- 
cesses. 

The points in Fig. [9] (bottom) concentrate to the coordinate 
origin, showing that in many quadruples the inner mass ratios and 
inner periods are similar, as in the e Lyr. The top plot of Fig. [Tol 
shows that the period disparity between the inner sub-systems in 
quadruple stars is less when Pl is short. This is expected because 
inner periods are comprised between few hours (contact binaries) 
and the dynamical stability limit, — 0.7 < log Ps < log Pl — 0.7. 
At shorter Pl this range is smaller, hence the inner periods are 
forced to be more similar, | log(Ps,i/Ps,2)| < log Pl ■ This range, 
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shown in Fig. [10] (top) by the dotted lines, is wider than the dis- 
tribution of the points by roughly 2 dex. Therefore the similarity of 
the inner periods is genuine and not caused by the trivial dynamical 
constraints. 

The numbers of quadruple systems with positive and negative 
| log(Ps,i/Ps,2)| are statistically similar, 38 and 43 respectively, 
so the periods of the more and less massive inner sub-systems do 
not differ systematically. The migration process that determined 
these periods should not depend on the mass. The inner periods 
are similar to within 1 dex in 42 systems, i.e. in 52 ± 8% of the 
quadruples. The relation between the ratio of the inner periods and 
the outer mass ratio is further explored in Fig. [10] (bottom). As a 
rule, sub-systems with similar inner periods also have similar total 
masses (larger qLi)- There are 34 quadruples (42% of the sample) 
with | log(Pgi/Pg2| < 1 and qn > 0.5, i.e. resembling e Lyr. 
Conversely, in the systems with dissimilar masses (low qn), the 
ratio of inner periods \og(Ps,i/ Ps,i) can take large positive or 
negative values with equal probability. 



6 SUMMARY AND DISCUSSION 

The properties of catalogued triple and quadruple systems are: 

(i) The distributions of periods in the inner sub-systems of 
quadruple and triple stars are similar and bimodal. They are dif- 
ferent from the period distribution of pure binaries. There is no 
correlation between inner mass ratios and inner periods. 

(ii) The distributions of periods and mass ratios in the outer 
sub-systems of quadruple and triple stars are different. The me- 
dian outer mass ratio in triples is 0.39. It does not depend on the 
outer period, which has a smooth distribution. In contrast, the outer 
periods of 25% quadruples concentrate in the narrow range from 
10 yr to 100 yr, all these tight quadruples have qn > 0.6 (possibly 
a selection effect), while the periods of their two inner binaries are 
similar. 

(iii) Quadruple systems of e Lyr type with comparable inner pe- 
riods and component masses are common. In 42% of the quadru- 
ples the outer mass ratios are above 0.5 and the inner periods differ 
by less than 1 dex. 

(iv) The outer and inner mass ratios in triple and quadruple stars 
are not mutually correlated. In 13% of quadruples both inner mass 
ratios are above 0.85 (double twins). 

(v) The values of the inner and outer orbital angular momenta 
(or periods) in wide (Ps > 30 d) triple and quadruple systems show 
some correlation, the ratio of outer-to-inner periods is mostly com- 
prised between 5 and 10 4 . In the systems with small period ratios 
(low hierarchy) the directions of the orbital spins are correlated, 
while in the systems with large ratios they are not. 

Statistical properties of multiple stars lead to several conclu- 
sions regarding their formation mechanisms. 

1. N-body dynamics is not the dominant process of 
multiple-star formation. Triple stars produced in N-body de- 
cay (with or without accretion) have distinctive properties: a nar- 
row period distribution (width of 1-2 dex), moderate period ratios 
Pl/Ps ~ 10 (not too far from the dynamical stability limit), high 
eccentricitie s of outer orbits distributed as /(e) = 2e or steeper 
[cf. Fig. 6 in Delgado-Dona te. Clarke. & Bate! ( 120031) and Fig. 6 in 
ISterzik & Tokovininl j2002)l, and low outer mass ratios qL3 < 0.2. 
Quadruple stars with 2+2 hierarchy are produced by this process 
only exceptionally. In contrast, in real multiple stars the p eriod ra- 
tios can be large, the outer eccentricities are moderate |shatsky 



l200ll : lTokovrnir]|2004 see also Table lATl . and the outer compan- 
ions are often massive (81% have qL3 > 0.2). Dynamical disrup- 
tion of small-N cluster s leaves too many single s tars and too few 
binaries and multiples dGoodwin & Kro upa 2005). 

The existence of outer weakly bound components or sub- 
systems, as in the cases of a Gem and 88 Tau (at projected dis- 
tances 1145 and 3500 AU, respectively), speaks against their dy- 
namical origin. If these multiples were parts of larger clusters, the 
stellar density must be low enough for the survival of the outer 
components, hence much too low to explain the dynamical forma- 
tion of the close inner sub-systems where larger binding energies 
are involved. 

2. The rotationally-driven (cascade) fragmentation sce- 
nario is promising. It naturally produces quadruple stars of 2+2 
hierarchy. As the motions in the pre-stellar cores are predominantly 
large-scale, the net angular momentum is sufficiently high, pre- 
venting disruptive N-body interactions between sub-systems. The 
net angular momentum is conserved in the orbital motion of the 
outer sub-system. The eccentricity of the outer orbit in this scenario 
should be moderate, the system is hierarchical and dynamically sta- 
ble. 

Suppose that two stellar pairs are formed from the fragments 
of one core. The initial separation between the components in each 
pair must be of the order of 10 2 to 10 4 AU. The proto-binaries inter- 
act with the surrounding gas, create spiral waves and transfer the or- 
bital angular momentum outwards, be coming tighter or even merg- 
ing in to single stars (disk migration, iKrumholz. Klein. & McKeg 
2007). The merger must not happen in the case of quadruples, how- 
ever. Each sub-system evolves independently of the other system 
until its closest reservoir of gas is exhausted. As the directions of 
the net angular momentum in each sub-condensation can differ, the 
orbital spins of the sub-systems are not necessarily co-aligned with 
each other. 

Accretion of the more distant envelope proceeds on a longer 
time scale, when the sub-systems already "fall" onto each other 
and end up forming a quadruple. The matter from the outer en- 
velope settles in a circum-quadruple disk. The result resembles a 
quadruple system GG Tau, where the circum-quadruple disk slowly 
feeds the circum-binary disk around the most massive sub-system 
Aab. Some components may form as fragments in accretion disks 
around other, more massive stars. The delayed formation of such 
components means that less mass is available for accretion, possi- 
bly creating binaries with low q. The AB quadruple in a Gem could 
be produced by such a process. 

3. Migration shortens both inner and outer orbits0 It is 
likely that there are several different types of migration. Accretion 
onto a binary and associated braking by the massive circum-binary 
disk is one such process confirm ed by hydrodynamical simulations 
( iBate. Bonnell, & Bromrr] r2002). Accretion usually increases the 
mass ratio and can produce twins with q ~ 1. However, we find 
no correlation between mass ratios and periods and infer that the 
dominant migration mechanism producing close binaries is not ac- 
companied by the accretion of substantial mass. An interaction with 
circum-binary disk is a promising candidate for such alternative 
migration mechanism. Such interaction also involves accretion, but 
possibly at a much reduced level not affecting the mass ratio. It is 
not clear whether migration can lead to a merger. If this is possi- 



2 As pointed out by the Referee, fragmentation at high densitie s may f orm 
close binaries or multiples without migration, cf. Machida et al] bOOijl) . 
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ble, massive tertiary components in triple stars could be mergers in 
former 2+2 quadruples. 

The gap in the distribution of the inner periods can be tenta- 
tively related to two distinct migration mechanisms. Presumably, 
accretion-induced migration acts at large separations, shortening 
the inner (and sometimes outer) periods to ~ 30 y (separation 
~ 10 AU) and creating a broad maximum in the period distribu- 
tion at log Pg ~ 5. The second migration mechanism shortens in- 
ner periods even further (mostly below 30 d) and produces the gap, 
but it does not affect the mass ratio. The second migration must be 
less efficient in pure binaries, so their perio ds are longer compare d 
to the inner sub-systems of multiple stars dTokovmin e t al. 2006). 
This explanation of the period gap is still highly hypothetical. 

In the case of tight quadruples, a significant shrinking of 
both inner and outer orbits must have happened. Most likely, 
the outer orbit shrinks by accretion from the envelope, increas- 
ing the outer mass ratio. This could explain high outer mass ra- 
tios of tight quadrupl e s (Fig . |7J and matches the simulations of 
iDelgado-Donate et alj d2004D . The gas accreted by the inner bi- 
naries inherits its angular momentum from the outer orbit, so the 
similarity of the inner periods in tight quadruples is a natural con- 
sequence of such process. We also expect the alignment of inner 
and outer angular momentum vectors, which is the case in some, 
but not all triple and quadruple stars. Accretion of even a small 
amount of gas can efficiently brake the binary if the angular mo- 
mentum of the gas is grossly mis-aligned with respect to the binary, 
as must happen in quadruple and triple stars with mis-aligned (e.g. 
counter-rotating) inner and outer systems. A tight quadruple with 
mis-aligned inner orbits, such as 88 Tau, could result from this pro- 
cess. It is noteworthy that both triples and quadruples with large 
Pl I Ps ratios show a tendency to mis-aligned orbital spins. 

To preserve dynamical stability, migration in the in- 
ner sub-systems must happen before or proceed faster than 
in the outer orbits. In the hydrodynamical simulations of 
iBate. Bonnell. & Bromml d2002l) some multiple systems become 
dynamically unstable because their outer orbits shrink too fast due 
to accretion and disk braking. We may be witnessing such situation 
in GG Tau which interacts wit h its circum-quadruple disk and is 
now close to the stability limit dBeust & Dutrevl l2006). A fraction 
of present-day triple stars may have been produced by dynamical 
decay of 2+2 quadruples. The triple star A Tau with Pl = 33 d 
and Pl / Ps = 8 could result from the outer-orbit migration which 
stopped just before the dynamical breakup. 

4. Kozai cycles with tidal friction (KCTF) is not the domi- 
nant process in producing close sub-systems. For long initial Ps, 
this mechanism "turns on" in a very narrow range of relative incli- 
nations near 90°. Moreover, at large Pl/Ps ratios as commonly 
found in triple systems the Kozai effect becomes too weak and is 
perturbed by other factors such as relativistic precession. Hence, a 
population of inner systems with sufficiently short Ps must be cre- 
ated by some other proce ss before the KCTF can furt her shorten 
some periods to few days (Fabrvckv & Tremaine 2007). This said, 
inner binaries with periods of ~ 1 d must be produced late, when 
the components have already settled on the MS, either by KCTF or 
by magnetic braking. 

Figure QT] illustrates the scenario of quadruple-star formation 
with &Pl—Ps diagram where some individual systems are located. 
The area of long inner and outer periods can be identified with the 
systems that did not evolve significantly after fragmentation. Oth- 
erwise, the points migrate to the lower left in this plot, faster in Ps 
than in Pl to stay always above the dynamical stability limit. 
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Figure 11. Schematic evolution of quadruple stars. The grey-scale back- 
ground shows the smoothed version of the Pl — Ps plot for quadruple 
systems (Fig. [5). 

Cascade fragmentation with subsequent orbit evolution is a 
promising explanation for the origin of multiple stars, including 
tight quadruples. This scenario remains very sketchy, however, as 
many details are still not clear, many alternatives remain to be ex- 
plored. From the theoretical side, simulations of accretion onto a 
multiple system with non-coplanar orbits are needed. Observation- 
ally, establishing the unbiased statistics of triple and quadruple stars 
in well-defined samples, such as nearby dwarfs, will be a crucial 
step in understanding the formation on these objects. We hope that 
the present catalogue-based study will stimulate such projects. 
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APPENDIX A: TIGHT QUADRUPLES 

A group of tight quadruples with log Pl < 4.5 is distinguished in 
the period-period plot (Fig. [3}- They also tend to have comparable 
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Table Al. Tight quadruples 
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masses and periods of the inner sub-systems. These properties are 
already obvious at log Pl < 5.4, so the exact upper limit Pl to 
distinguish tight quadruples from the rest is not known. We set this 
limit at Pl = 100 y (log Pl < 4.56) and study the properties of 
such quadruples in more detail. 

Table lAll lists all 20 known quadruples with Pl < 100 y: the 
WDS codes, alternate identifications, outer periods, apparent semi- 
major axes (or separations) of the outer orbits at, their eccentrici- 
ties ei, and primary masses A/i ,i . Only 3 out of 20 quadruples still 
have no computed outer orbits (missing eccentricities in Table lAlt . 
In 7 systems there are 5 or 6 known components, like in 88 Tau. 
This information is complementary to the main Table 3. 

The separations in the outer orbits are sub-arcsecond in 18 
systems, so the study of the inner orbits did not rely on spatially 
resolving the outer pairs, but rather on the spectroscopy of com- 
ponents in the combined light. Short outer periods can be easily 
discovered by spectroscopy, therefore the sharp drop in the number 
of quadruples at logP^ < 3.5 must be genuine, not a selection 
effect. T. Pribulla (private communication) confirms that in a large 
sample of triple stars discovered spectroscopically the outer peri- 
ods are mostly long and the tertiary companions are mostly single, 
so the quadruple VW LMi with Pl = 355 d stands out as a unique 
system. 



APPENDIX B: MASS RATIOS WITH INDEPENDENT 
PAIRING 

Is is established that the mass ratios in binary stars do not 
correspond to the random selection of component masse s from 
the initial mass function (IMF) (e.g. Halb wachs et al] 120031 : 
iKouwenhoven et al]|2007l) , the component masses rather tend to 
be comparable. The same tendency is expected in the multiple sys- 
tems. Nevertheless, it is instructive to simulate mass ratio distribu- 
tions resulting from the independent pairing and to compare th em 
with the data. We used the numerical recipe of lFished d2004h to 
generate the masses. Only triples with Mi > M m i n and quadruples 
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Figure Bl. Mass ratios in 865 simulated triples (top) and 463 simulated 
quadruples (bottom) with independent component masses selected from the 
IMF. Compare with Fig. [5] 



with Mi,i > Af m in, with Mmin = 0.8Mq, were retained from a 
larger pool of simulated systems dominated by small masses. This 
is done to mimic our samples, although the distribution of simu- 
lated primary masses is more strongly concentrated towards M m i n 
than the actual masses (Fig. |2). The plots of the simulated mass 
ratios are shown in Fig. lBll The difference with the actual distribu- 
tions is too strong to be attributed solely to selection. If the masses 
in the real multiple stars were indeed independent and we were 
missing a huge number of multiples with low mass ratios, it would 
mean that the actual frequency of multiple systems is many times 
higher than observed now. Systems with similar components are 
rare in the simulated multiples. The median simulated mass ratios 
are qLs, = 0.13, §L4 = 0.29. A loose correlation between the small 
values of qs and qLz seen in the Figure is explained by the thresh- 
old imposed on the primary mass, while the other components most 
frequently have similarly low masses near the IMF peak. 

This paper has been typeset from a TpX/ KTpX file prepared by the 
author. 



